YAMABE FLOW ON MANIFOLDS WITH EDGES 

ERIC BAHUAUD AND BORIS VERTMAN 

Abstract. This article studies local existence and uniqueness of Yamabe flow within a class 
of compact Riemannian manifolds with incomplete edge singularities. Our main analytic step 
is to establish parabolic Schauder-type estimates for the heat operator on a class of edge 
singular spaces. We apply these estimates to obtain local existence for certain quasilinear 
^^ . equations, including the Yamabe flow. 
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> ■ 1. Introduction and statement of the main result 

O 

On a compact Riemannian manifold, the Yamabe problem asserts that every conformal 

class of metrics contains a representative of constant scalar curvature. There are now several 

proofs of this fact. The first proof, commenced in [Yam60] and continued by [Tru68, Aub76, 

Sch84] used the calculus of variations and elliptic partial differential equations; see [LP87] 

for a survey. Another proof uses the geometric Yamabe flow: 

d t g = -scal(g(t)) ■ g, 



CO 

in 



(^ #(0) = £init, 

(and appropriate normalizations) to converge to constant scalar curvature metrics; see 
[Brell] for a recent overview. 

It is natural to wonder to what extent the Yamabe problem holds in other settings. There 
has been recent work in understanding this problem on singular manifolds with conic and 
more general incomplete edge metrics. See the work of Akutagawa-Botvinnik [AB03] for the 
conic case, cf. also Jeffres-Rowlett [JR10]. Further, see Akutagawa-Carron-Mazzeo [ACM11] 
for the Yamabe problem on very general stratified spaces. Note that this is not what is 
commonly called the singular Yamabe problem which asks that given an appropriate closed 
subset r C M find a complete metric of constant scalar curvature on M\T. All of the papers 
mentioned above attack the problem from an elliptic PDE point of view. 

Ricci flow on surfaces starting at singular metrics has been studied by various authors, see 
the recent review by Isenberg-Mazzeo-Sesum [IMS 11]. One expects uniqueness of the flow 
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to fail without specifying further boundary restrictions. Not only is it possible to obtain a 
Ricci flow that remains singular (see Mazzeo-Rubinstein-Sesum [MRS11] and Yin [YinlO]), 
but one may also obtain a solution to the flow starting at a singular metric that becomes 
instantaneously complete (see Giesen- Topping [GTlOb, GTlOa]) or smooths out the cone 
angle (see Simon [Sim02]). 

In this paper and its successor we are interested in the Yamabe flow on spaces with in- 
complete edge singularities that preserves the singular structure, as a first step to study the 
Yamabe problem from the geometric flow perspective within the setup of singular spaces. 
Our work is for general dimension n > 2 and allows for more general incomplete edges. 

On closed manifolds local existence of Yamabe flow is a basic consequence of the classical 
existence theory for parabolic partial differential equations, see Ladyzhenskaya-Solonnikov- 
Uraltseva [LSU67]. In the presence of conical or edge singularities the analysis is complicated 
by the fact that both the operators and the scalar curvature are singular. Thus geometric 
and analytic conditions must be imposed to even make sense of the flow. We mention that 
in the setup of isolated conical singularities, existence and regularity of solutions to the 
inhomogeneous heat equation has recently been addressed in a recent preprint of Behrndt 
[Behl 1]. Another approach to estimates for conical singularities, and different from the one 
we give here, is given in the forthcoming paper by Mazzeo-Rubinstein-Sesum [MRS11]. 

Our main result is as follows; see Theorem 1.8 below for the precise statement specifying 
the class of metrics we consider. 

Theorem 1.1. There exists a solution to the Yamabe flow starting within a class of com- 
pact Riemannian spaces with admissible simple edge singularities that remains asymptotically 
admissible for a short-time. 

To give a concrete example of a metric for which our theorem applies, one may take a 
direct product of a cone over S>^ with the round metric and any compact manifold. Our 
theorem also applies to any sufficiently high order perturbation of this metric and certain 
non-product cone bundles. 

The remainder of the introduction is organized as follows. We next discuss the class of 
incomplete edge metrics. We then describe our parabolic Holder spaces and state our main 
analytic result that provides an analogue of classical parabolic Schauder theory in the singular 
setup. Finally we state our short-time existence result for the Yamabe flow precisely. 

1.1. Simple edge spaces. Let M be a compact stratified pseudomanifold, with top- 
dimensional open and dense stratum M and a single lower-dimensional stratum B, which is 
a compact smooth manifold, by the axioms on the differential and topological structure of 
stratified spaces. 

Moreover, B has a neighbourhood <fy C M, a radial function x : % — > M + and a smooth 
bundle projection <ft : % — y B, which is a submersion on ty = % D M, such that the 
preimages <p~ l {q) H % are all diffeomorphic to open cones (0, 1] x F over a compact smooth 
manifold F, and where the restriction of x to each fibre 4>~ 1 {q) is a radial function on that 
cone. The level set Y := x (1) is the total space of a fibration : Y — y B with fibre F, 
smooth and compact, and corresponds to the regular boundary of the neighbourhood % . 

We fix dimensions for the remainder of the paper. We let m = dim M, b = dim B and 
/ = dim F, note that 

m=l + f + b. 
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We say that a metric go on M is a rigid incomplete edge metric if it is smooth and 
Riemannian away from fy , and has the following form in ^: there is a smooth Riemannian 
metric f) on B and a symmetric 2-tensor k on Y, restricting to Riemannian metrics on fibres 
F, such that 

g Q \w = dx 2 + (j)*l) + x 2 k. 

Locally in trivializing neighbourhood of the fibration <p, this describes ty simply as the 
cone bundle over B with open truncated cones (0, 1] x F as fibres, with the metric induced 
by go on eac h fibre cj)~ 1 (q) being an exact warped-product conic metric. 

More generally, we call g an incomplete edge metric if g = g + h where go is rigid in 
the sense above and where h is a smooth symmetric 2-tensor such that \h\ go is smooth on 
[0, l) x x Y and vanishes at x — 0. We refer to any compact stratified Riemannian space with 
the structure just described, and with an incomplete edge metric, as a simple edge space 
(M, g) with edge data (B, F, <p). 

Among the rigid incomplete edge metrics there is a slightly more restricted class of admis- 
sible edge metrics go, where <fi : (Y, g\ Y ) — > (B, f)) is a Riemannian submersion with respect 
to g . Recall, if p G Y, then T P Y splits into vertical and horizontal subspaces, TYY © T^Y , 
where by definition T^Y is the tangent space to the fibre of <fr through p and T^Y is the or- 
thogonal complement of this subspace. The submersion condition implies that the restriction 
of the tensor k to T^Y vanishes. 

We now present the class of metrics on which we can establish Schauder-type estimates 
for the Laplacian. 

Definition 1.2. Let (M,g) be a simple edge space with edge data (B,F,<p) and a singular 
neighbourhood ^cMof the edge B, such that the incomplete edge metric g is of the form 
9 — 9o + h with 

9o\w — dx 2 + (j)*t) + x 2 k. 
We call (M, g) a feasible edge space, if the following five conditions are satisfied 

(i) \h\ go vanishes to second order at x — 0, i.e. \h\ go = 0(x 2 ) as x — > 0, 

(ii) is a Riemannian submersion with respect to go, 
(iii) the restrictions (F, «|^-i(n,)), b G B are isospectral, 
(iv) the Laplacian A go preserves the space of sections that are fibrewise constant over °i/ , 

(v) the Laplacians A K ^ associated to (F, «|^-i(6)) for any b G B, satisfy 

SpecA Kife \{0} > dimF. 

We refer to the simple edge spaces satisfying the first two conditions of Definition 1.2 
as asymptotically admissible. This class of singular spaces has already been considered in 
[MV11]. The first two assumptions comprise the basic geometric and analytic conditions 
required for analysis of edge spaces. The more interesting geometric restriction is given by 
the last three conditions in Definition 1.2. We discuss these in the examples below. 

Example 1.3. Let (M,g) be an asymptotically admissible edge space with edge data 
(B,F,<f)). 

(i) The conditions (iii) and (iv) in Definition 1.2 are satisfied in case of a trivial fibration 
_1 ({1}) = B x F. This in particular holds in the special case where B is a point 
and the edge reduces to an isolated conical singularity. 

(ii) Given a Riemannian submersion <fi as above, a vector field H on Y is projectable 
if there exists a vector field H B on B where <fi*H = H B . In [Mul09, Satz 2.9], 
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Jorn Miiller proved a sufficient condition for Definition 1.2 (iv) is that the mean 
curvature H of the fibres F of _1 ({1}) is projectable. Two sufficient conditions for 
projectability are that the fibres F to be totally geodesic ([Mul09, Section 2.12]), 
in which case fibres F are isometric and further Definition 1.2 (iii) is also satisfied, 
see [Bes08, Theorem 9.56]. A second condition is that (M,go) is locally a warped 
product, see [Bes08, Theorem 9.104] and [FIP04]. 

The conditions (iii) and (iv) in Definition 1.2 provide a restriction of the geometry of the 
fibration 0, whereas the final condition (v) in Definition 1.2 on the spectrum of fibres is 
analytic and discussed in the examples below. 

Example 1.4. Let (M m ,g) be an asymptotically admissible edge space with edge data 
(B b , F* , <fi) and a singular neighbourhood % C M of the edge B. We make the dependence 
of g on the metric structures on fibres F and the base B explicit by writing g = g(t),K). 
Assume that the Laplacians A K ^ associated to (F, «|<£-i(&)) for any b G B are isospectral and 
A > is the lowest non-zero eigenvalue. 

(i) For any c > y/f/Xo the edge space (M, g(t), c~ 2 k)) satisfies Definition 1.2 (v). 
(ii) In case of dimM = 2 and F = S , the condition (v) of Definition 1.2 is satisfied 

whenever the singular neighbourhood f cMof the edge B is a fibration of cones 

of a fixed cone angle strictly less than 7r/4. 

1.2. Mapping properties of the heat operator on Holder spaces. The proof of short- 
time existence for Yamabe flow is based on a careful analysis of the mapping properties of 
the heat operator acting with respect to Holder spaces, defined with respect to an incomplete 
edge metric. Consider a feasible edge space (M , g) and recall the notation fixed above. 

Definition 1.5. The Holder space A a 'a(M x [0,T]),a G (0,1), is defined as the space of 
functions u(p,t) that are continuous on M x [0,T], such that their a-th Holder norm 

/ \u(p,t)-u(p',t)\ \ 

\\u\\ a := Moo + sup — — -pr < oo, 

\d M (p,p') a + \t-t'\2 J 

where cIm(p,p') represents the distance between p,p' G M with respect to the incomplete 
edge metric g and in local coordinate charts in the singular edge neighbourhood % may be 
equivalently defined by 



d,M((x,y, z),{x',y', z')) = \/\x — x'\ 2 + (x + x') 2 \z — z'\ 2 + \y — y'\ 2 . 

Definition 1.6. Consider the Laplacian A g of (M, g) and the edge vector fields V e introduced 
in [Maz91] and reviewed at the beginning of Section 2. Let V' e be a finite set of local generators 
of V e and put 

f:=RA 9) i-V| VeV e }. 

Then the higher order Holder spaces are defined as follows 

A 2k+a,k+f( Mx j 0)T ]) = { MG A a 'f |Aju,AoA^nG A Q 't(M x [0,T]), 

Xe@,j = 0,...,k-1}, 

with differentiation is a priori in the distributional sense and the Holder norm 

fc-i fc-i 

||«||2*+a = \\ u \\a + ^2 II A 9 n H° + 5Z 5Z ^ X ° A 3 M H Q - 
j=0 j=0 X€$> 
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Our first result discusses the mapping properties of the heat operator on simple edge spaces, 
and essentially correspond to the classical parabolic Schauder estimates. 

Theorem 1.7. Let (M m ,g) be a feasible edge space with edge data (B b , F* ,(f>) and a singular 
neighbourhood % C M of the edge B. The incomplete edge metric is given by g = g + h 
with 

9o\<% = dx 2 + 4>*i) + x 2 k. 

Let \ Q > f be the smallest non-zero eigenvalue of A Kb ,b G B. Put 



(/-l) ,/ (/-l) 2 
«o := ^ h V 4 h A - 1 > 0. 

Denote by A g the Friedrichs extension of the Laplacian on (M,g). Then for a G 
(0, min{ao, 1}) the heat operator e~' As acts as a bounded convolution operator with respect to 
the Holder spaces (k G N ) 

e -tA g . A 2k+a,k+f (M x [Q)T]) ^ A 2(k+l)+a,(k+l)+f (M x ^^ 
e -tA g . A 2k+ Q ,k+f ^ x [ 0)T ]) -> v^A2k+a,k+f ( M x [ 0j T]), 

As elaborated in Example 1.4, the relation A > dimF can be achieved using a rescaling 
of the metric k, and corresponds to choosing small cone angles in the surface case. This 
restriction guarantees ao > 0. For dimM = 2 and the singular neighbourhood ^ cMof 
the edge B being a fibration of cones over F = S 1 of a fixed cone angle 6 < 7r/4, we have 
ao = cotan^ — 1 > 0. In the rest of the paper ao and a G (0, minjao, 1}) are fixed. 

We take this opportunity to mention a related paper. In joint work with Emily Dryden 
[BDV11], we have also obtained mapping properties of the heat operator for the homogeneous 
Cauchy problem for the heat equation using time-weighted (spatial) Holder spaces adapted 
to the incomplete metric. The estimates we obtained there are more naturally applied to 
semilinear parabolic problems. 

1.3. Statement of the main result. We now return to the Yamabe flow given in equation 
(1.1). Let (M m , <7i n it) be a feasible edge space. Straightforward computation shows that scalar 
curvature becomes unbounded near the edge x = 0, so further restrictions are required to 
ensure the flow exists. In addition to the feasibility hypothesis we will require 

scal( feit )GA 2fe+Q (M,3 init ), (1.2) 

for some k G N, where A 2fc+Q (M,# init ) C A 2k+a ' k +f (M x [0,T]) is the subset of all time- 
independent functions of A 2k+a ' k+ 2"(M x [0,T]). We discuss the hypothesis (1.2) for lowest 
regularity k = 1 more carefully in Section 5 but for now we point out that it implies that the 
obstruction scal(ft) = /(/ — 1) is satisfied, cf. [AB03]. 

Returning to the Yamabe flow, since the flow preserves the conformal class of g m it we write 

g(t) = e 2 <^ ■ flfa*, 



to 



which transforms the Yamabe flow to a scalar equation for u 

d t u = -V 2 " (2(m - l)A ginit u - ( m -2)(m- l)\Vu\ 2 + scal(^ nit )) , u(p,0) = 0, 

where A 9init is the Laplacian of (M, g; n it). This is a quasilinear equation for u and our main 
result concerns local existence and uniqueness of its solutions. 
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Theorem 1.8. Let (M m ,g init ) be a feasible edge space such that scal(g ini t) G A 2k+a (M, g init ) 
for some k G N. Then the (transformed) Yamabe flow equation admits a unique solution 
u G A 2k+Q > k +§(M x [0,T]) for some T > and a G (0,a ). 

Note that e 2u g- m n is again an asymptotically admissible simple edge metric after a suit- 
able transformation of the defining function. Thus, Theorem 1.8 indeed asserts short time 
existence of Yamabe flow within a class of compact Riemannian spaces with simple edge 
singularities. 

The outline of the remainder of this paper is as follows. In Section 2 we discuss the 
asymptotics of the heat kernel of the Laplacian. In Section 3 we derive the analogue of 
Schauder estimates. In Section 4 we prove a short-time existence result for certain quasilinear 
parabolic equations, and in Section 5 we condition the Yamabe flow to apply the result of 
Section 4. 

Due to the length of this paper we decided to postpone a discussion of long-time existence 
and convergence results for the Yamabe flow for these metrics. This will be done in a 
forthcoming manuscript. 
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authors would like to thank Simon Brendle for helpful remarks. The second author gratefully 
acknowledges financial support by the German Research Foundation DFG as well as the 
Hausdorff Institute in Bonn, and also thanks Stanford University for hospitality. 



2. Asymptotics of the heat kernel on simple edge spaces 

Let (M, g) be a simple edge space with edge data (B, F, 0) and a singular neighbourhood 
f C Mof the edge B. We pass from the singular space M to its resolution M, which is 
a manifold with boundary Y = dM, obtained by replacing % , which is a bundle of cones, 
with the associated bundle of cylinders, where each fibre is now [0, l) x x F. This resolution 
process is described in greater detail in [Maz91]. 

Consider local coordinates (x, y, z) on M near the boundary Y = dM, where x is the 
radial coordinate, y is the lift of a local coordinate system on B b and z restricts to a local 
coordinate system on each fibre F* . The class of edge vector fields V e on M are those which 
are smooth and tangent to the fibres of Y at at dM. In this local coordinate system, any 
element of V e can be written as a sum of smooth multiples of the basic generators xd x , xd yi 
and d~., which we write as 



V e = C°° -span {xd x ,xd yi ,...,xd Vb ,d Zl ,...,d Zf } . 

The edge vector fields define the class of differential edge operators Diff*(M). By definition, 
L G Diff*(M) if it can be written locally as a sum of products of elements of V e , with 
coefficients in C°°(M). Thus 

L= 5Z a ha,p{x,y,z){xd x ) 3 {xd y ) a d^ 

j+\a\ + \p\<n 
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with each a^ a fi smooth. More generally, if L acts between sections of two vector bundles, 
then L has this form with respect to suitable local trivializations, where each a,- jaj/ g is matrix- 
valued. The operator is called edge elliptic if its edge symbol 

e a n (L)(x,y,z;^r],C)-= J^ a j,aA x >V> z )t 3r l a t P 

j+\a\+\0\=n 

is nonvanishing (or invertible, if matrix- valued) , for (£, rj, () ^ (0, 0, 0). This has an invari- 
ant meaning as a function on the so-called edge cotangent bundle e T*M which is homogeneous 
of degree n on the fibres. The asymptotic structure of solutions to an elliptic differential edge 
operator L is encoded in its indicial operator, which acts on functions on R + x F, for any 
fixed point y G B 

H L )vo : = J2 a J,oA Q ^yo,z)(sd s yd^. 

j+\P\<m 

This is equivalent, by taking the Mellin transform in the s (g 1R + ) variable, to a family of 
holomorphic operators on F, 

k( L )yo= J2 a ifiA Q iyo,z)( J dz- 

j+\/3\<m 

Values of ( G C for which I^(L) yo is not invertible (acting on L 2 (F)) are called indicial roots 
of L; they are analogous to eigenvalues for this operator family on the compact manifold F. 
For the scalar Laplacian A of a feasible edge space (M, g), x 2 A G Diffg(M). The fibration 
<p : Y — > B is a Riemannian submersion with respect to the rigid part go of the edge metric 
9 — 9o + h, where go \ % — dx 2 + 0*f) + x 2 k and \h\ go = 0(x 2 ) as x — > 0. Consider the 
orthogonal splitting of TY into vertical and horizontal subbundles 

TY = kerd0©(kerd0) ± , 

where dip : (ker d<f)) ± — > B, is an isomorphism and the horizontal subbundle ker d<f) is 
annihilated by k. Under the identification (ker dip) 1 - = B the Laplacian A is given over the 
singular neighbourhood % by the following expression 

x 2 A \ <% = -{xd x f - (dimF - l)xd x + A F + x 2 A B + xV, V G V e , 

where A^ is the pullback of the Laplacian of (F, K\<p-in,)) to fibres of the fibration <p : Y — > 
B, and A^ is the lift of the Laplacian for (B, 1)). The higher order term V G V e is comprised 
of the contributions from h as well as the second fundamental form and the curvature of 
the fibration <p. The Mellin transform of the corresponding indicial operator is given by the 
holomorphic operator family on F 



/ c (x 2 A) m) = -C 2 -(dimF-l)C + A 



KjW 



(F, /■£ | </,-!(&)), b G B, are spectrally equivalent and hence indicial roots of x A are y - 
independent. 

Denote the Friedrichs extension of the Laplacian on (M, g) again by A. We denote its heat 
operator by e~ tA and the corresponding heat kernel by H. The heat operator of A acts as 
an integral convolution operator on f(t, •) G @(A),t > 0, 



-tA 



*f(t, P ) 



[ f H(t-t,p,p)f(t,p)dvol g (p)dt, (2.1) 

JO J M 
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and solves the inhomogeneous heat problem 

(d t + A)u{t,p) = f(t,p), u(0,p) = 0, 

for any /(£,•) G @(A),t > 0. The heat kernel can be viewed a priori as a distri- 
bution on Mf L = IR + x M 2 , with the local coordinates near the corner in M% given by 
(t, (x, y, z), (x, y, z)), where (x, y, z) and (x, y, z) are the coordinates on the two copies of M 
near the singularity. The kernel H(t,(x,y,z),(x,y,z)) has non-uniform behaviour at the 
submanifold 

A = {(t = 0, (0,y,z), (0,y,z)) eR + xdM 2 \y = y}, 

and the diagonal D, which requires an appropriate blowup of the heat space M 2 , such 
that the corresponding heat kernel lifts to a polyhomogeneous distribution in the sense of 
the following definition. 

Definition 2.1. Let 2U be a manifold with corners, with all boundary faces embedded, and 
{(Hi, Pi)}f = i an enumeration of its boundaries and the corresponding defining functions. For 
any multi-index b = (&i, . . . , bjy) G C N we write p b = p± . . . p^ . Denote by Vf,(2U) the space 
of smooth vector fields on 2U which lie tangent to all boundary faces. A distribution u on 2U 
is said to be conormal, if w G p 6 L°°(QU) for some b G C N and ]/]_.., VpU G p b L°° (Q2J) for all 
Vj G Vfc(2U) and for every £ > 0. An index set Ei = {(7,p)} C C x N satisfies the following 
hypotheses: 

(i) Re(7) accumulates only at plus infinity, 

(ii) For each 7 there is P 7 G No, such that (j,p) G Ei for p < P 1 < 00, 
(iii) If (j,p) G Ei, then (7 + j,p') G Ei for all j G N and < p' < p. 

An index family E = (Ex, . . . , E^) is an iV-tuple of index sets. Finally, we say that a conormal 
distribution w is polyhomogeneous on 2U with index family E, we write u G ,g^ g (2U), if w 
is conormal and if in addition, near each Hi, 

(7>P)6Ei 

with coefficients a 7P conormal on i/j, polyhomogeneous with index Ej at any if, D Hj. 

The basic underlying idea of the heat-space blowup is to capture the well-known scaling 
property of the cone and the resulting polyhomogeneous behaviour of the heat kernel as 
its entries approach certain submanifolds of M% from the various angles. This is done by 
"blowing" up these submanifolds, parabolically in time-direction, a procedure introduced by 
Melrose, see [Mel93]. 

The blowup procedure corresponds to introducing certain (polar) coordinates on M\, such 
that the heat kernel admits polyhomogeneous asymptotic expansions at A and D in these 
coordinates, together with a unique minimal differential structure with respect to which these 
coordinates are smooth. To get the correct blowup of M% we first parabolically blow up the 
submanifold 



A = {(t, (0,2/,*), {0,y,z)) eR + xdM 2 :t = 0,y = y}c M, 



2 
h- 



The resulting heat-space [M%, A] is defined as the union of M%\A with the interior spher- 
ical normal bundle of A in M\. The blowup [M%, A] is endowed with the unique minimal 
differential structure with respect to which smooth functions in the interior of M| and polar 
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coordinates on M\ around A are smooth. This blowup introduces four new boundary hy- 
persurfaces; we denote these by ff (the front face), rf (the right face), If (the left face) and tf 
(the temporal face). 

The actual heat-space blowup ^M 2 is obtained by a parabolic blowup of [M%, A] along the 
diagonal 

D ■■= {t = 0, x = x, z = z, y = y] C M 2 hl 

lifted to a submanifold of [M%, A]. The resulting blowup ^% is defined as before by cutting 
out the submanifold and replacing it with its spherical normal bundle. It is a manifold with 
boundaries and corners, visualized in Figure 2. 




Figure 1. Heat-space Blowup Jt^. 



The projective coordinates on ^^ are then given as follows. Near the top corner of the 



front face ff, the projective coordinates are given by 



P 



Vt, £ 



X 



p 



z 



X 



p 



u 



y-y 

■ 
p 



z, y, z, 



(2.2) 



where in these coordinates p, £, £ are the defining functions of the boundary faces ff, If and 
rf respectively. For the bottom corner of the front face near the right hand side projective 
coordinates are given by 



r 



t 



x' 



X 



X 



II 



y-y 



X 



z, x, y, z, 



(2.3) 



where in these coordinates r, s, x are the defining functions of tf, rf and ff respectively. 
For the bottom corner of the front face near the left hand side projective coordinates are 
obtained by interchanging the roles of x and x. Projective coordinates on ^£ near temporal 
diagonal are given by 



V 



Vt 



X 



S 



[x — X) 



Vi 



U 



y-y 

Vt ■ 



X[Z 



2) 



Vt 



x, y, z. 



(2.4) 



In these coordinates tf is the face in the limit \(S, U,Z)\—too,& and td are defined by 
x, 7], respectively The blow-down map (3 : e/#^ — > M\ is in local coordinates simply the 
coordinate change back to (t, (x,y,z), (x,y,z)). The heat kernel lifts to a polyhomogeneous 
conormal distribution on M 2 with product type expansions at the corners of the heat space, 
and with coefficients depending smoothly on the tangential variables. More precisely we 
have for asymptotically admissible edge metrics, notion introduced in [MV11], which are 
incomplete edge metrics, satisfying conditions (i) and (ii) of Definition 1.2. 
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Theorem 2.2. ([MV11]) Let (M m ,g) be a simple edge space with an asymptotically ad- 
missible edge metric g. Let the heat kernel associated to the Friedrichs extension of the 
corresponding Laplacian be denoted by H . The lift f3*H is a polyhomogeneous conormal dis- 
tribution on ^^ of leading order (— dimM) at the front face ff, leading order (— dimM) at 
the temporal diagonal td, and continuous at the left and right boundary faces. (3*H vanishes 
to infinite order at the temporal face tf and does not admit logarithmic terms in its asymptotic 
expansion at ff and td. 

The asymptotic behaviour of the heat kernel is in fact derived in [MV11] under a unitary 
transformation similar to [BS88], which changes the leading order behaviour at ff. The precise 
argument is outlined in [BDV11]. 

In fact, feasibility of the edge metric guarantees a finer result on the asymptotic behaviour 
of H at the left and right boundary faces of ^1- 

Proposition 2.3. Let (M m ,g) be a feasible edge space with edge data (B b , F* , <fi) . Then the 
lift (3*H is a polyhomogeneous distribution on ^£ with the index set at rf given in terms of 



E = {7 > | 7 = ~^— ^ + J (/ 4 1)2 + A, A G Spec A F }. 

More precisely, /3*H separates into two components /3*H' and f3*H" of leading order be- 
haviour (— m) and (— m + 1) at the front face, respectively; and if s is a defining function of 
the right face of the form (2.2) or (2.3) and E* = E\{0}, then the expansion of f3*H' and 
/3*H" takes the following form 

f3*H' ~ J2 s 7 a 7 (/3*^), s -> 0, 

7S-E 
oo oo 

f3*H" ~ J2 sla i(P*H) +J2J2 s 1+l a lt i{f5*H), s -+ 0, 

1=2 7G£* 1=1 

where ao(/3*H), 0,2(0* H) are constant in z. 

Proof. Recall the heat kernel construction in [MV11]. The initial approximate parametrix for 
the solution operator of C = (d t + A) is constructed in [MV11] by solving the heat equation 
to first order at the front face ff of ^^- The restriction of the lift (3*(tC) to ff is called 
the normal operator Ns{t£) at the front face and is given in projective coordinates (2.3) 
explicitly as follows 

N s (tC) = r(d T - dl - fs~ l d s + S - 2 A F , Z + Af ) =: r(d T + Af <*> + Af ). 

Nff(tC) does not involve derivatives with respect to (x,y,z) and hence acts tangentially to 
the fibres of the front face. Searching for an initial parametrix H , we solve the heat equation 
to first order at the front face, and note that 

N s (tC o H ) = N s (t£) o N S (H ) = 

is the heat equation on the model edge ^(F) x M. b w , where ^(F) = IR+ x F z denotes the 
model cone. Consequently, the initial parametrix H is defined by choosing Ns(H ) to equal 
the fundamental solution for the heat operator Ng(t£), and extending Nff(H ) trivially to a 
neighbourhood of the front face. Using the projective coordinates (r, s, u, z, x, y, z) near ff , 
see (2.3), we have 

N S (H ) :=H^ f \t,s,z,1,I)H r \t,u,0), (2.5) 
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where H R denotes the euclidean heat kernel on K 6 , and H^( F > is the heat kernel for the 
exact cone ^(F), as studied by [Che83], [Les97] and [Moo99]. Consequently, the index set 
Eq for the asymptotic behaviour of Hq at the right and left boundary faces is given in terms 



V(F) 



on the exact cone, i.e. 



of the indicial roots 7 for the Laplacian A s 

-y&E 

with the leading coefficient ao(H Q ) being harmonic on fibres and hence constant in z. Note 
that by condition (v) of Definition 1.2, any 7 7^ is automatically 7 > 1 + «o- The error of 
the initial parametrix Hq is given by 



P*(t£)H = [P*(tA g ] 
Due to g = go + h with \h\ an = 0(x 2 ) we have 



#0 =: Po. 



I.yo 



A„ - A, 



, fl ^ go - V X 6 V e 

Consequently in coordinates (r, s, w, z, x, y, z) near ff fl rf 



P = T3?/3*(Vi)H + f/3*(£A go ) - rAf ^ - rA 



F n 



As for the first summand, note that the asymptotic expansion of Hq as s — > starts with 
s°ao with no s°(logs) fc , k > 1 terms. The leading term s°ao is annihilated by the edge vector 
fields sd s and <9 Z , since a is constant in z, and its order is raised by sd u . Hence 

00 00 



1=1 



7 G£* /=o 



where the coefficients are of leading order (— m + 2) or higher in x at the front face. As for 



the second summand, note that 



p*(tA 



.</(} / 



rA 



V(F) 



rA? 



is of higher order in s and is of 



higher order in x, since its normal operator at the front face is zero by construction. We find 

"/r(£A,J-rA5 F) -rA? 



J 90/ 



S l Ci 



EE s7 ~ 1+ ^> ^°- 



1=0 7S-B* 1=0 

where the condition (iv) of Definition 1.2 implies that Cq is constant in z, and all coefficients 
are of leading order (— m + 1) or higher at the front face. Therefore Po is of leading order 
(— m + 1) at the front face and its expansion at rf is given by 

00 00 

i=0 7GE* Z=0 

where ao(Po) is constant in z. The next step in the construction of the heat kernel involves 
adding a kernel H Q to P^o, such that the new error term is vanishing to infinite order at rf. 
In order to eliminate the term s 7 a 7 in the asymptotic expansion of Po at rf, we only need to 
solve 

{-d 2 s - fs- l d s + s- 2 A F , z )u = s^r- 1 ^). (2.6) 

This is because all other terms in the expansion of tC at rf lower the exponent in s by 
at most one, while the indicial part lowers the exponent by two. The variables (r, u, x, y, 1) 
enter the equation only as parameters. The equation on the model cone is solved via the 
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Mellin transform and the solution is polyhomogeneous in all variables, including parameters 
and is of leading order (7 + 2). Consequently, the correcting kernel H' Q must be of leading 
order 2 at rf and of leading order (— m + 1) at ff, since Pq is of order (— m + 1) at ff and the 
defining function x of the front face enters (2.6) only as a parameter. 
In case of 7 = 0, the indicial equation (2.6) reduces to 

(-d 2 s -fs- 1 d s )u = s°(T- 1 a ), 

with a and the solution u = — s 2 a (2 + 2/) _1 both independent of z. Hence, the leading 
order term in the expansion of H' Q at rf is constant in z. Put Hi := H + H' Q , H' q is of leading 
order (— m + 1) at ff and expands near rf as follows 

00 00 

H' -JAK) + J2 5> 7+Zg ^o), s -► 0, 

1=2 *y€E* 1=1 

where we have explained that 02(^0) is constant in z. The error of our new heat parametrix 
Hi is now vanishing to infinite order at rf and is of leading order (— m + 1) at ff. 

In the following correction steps the exact heat kernel is obtained from Hi first by improving 
the error near td and then by an iterative correction procedure, adding terms of the form 
Hi o (Pi) k , where Pi := tCHi is vanishing to infinite order at rf and td. 

Extend s 7 a 7 (/7i) to a polyhomogeneous distribution on ^#^ with same index set at ff and 
If as Hi, vanishing to infinite order at tf and td, and a single term s 7 a 7 (iJi) in its expansion 
at rf. Then s 7 a 7 (ifi) o (Pi) k also makes sense, and since (3*(xd x — 7)s 7 a 7 (i/i) vanishes to 
infinite order at rf, so 

(3*(xd x - 7 )[s 7 a 7 (#i) ° ( p i) k ) ~ 0(s°°), s -> 0. 

Consequently s 1 a 1 {Hi) o (Pi) k has a single s 7 term in its expansion at the right face. 
Similarly, applying d z , we show that the coefficients in the expansion of s°ao(Hi) o {Pi) k and 
s 2 a 2 {Hi) o (Pi) k at the rf are constant in z. We deduce that (3*H = (3*H' + f3*H", with f3*H' 
of leading order (— m) at ff, (3*H" of leading order (— m + 1) at ff, and 

f3*H f ~ J2 ^a-iiPH), s -> 0, 

00 00 

f3*H" ~ J2 sla i^* H ) +Y.Y1 ^ +l a ltl (/3*H), s -)• 0, 

1=2 7S-B* 1=1 

where ao(fl*H), 0,2(0* H) are constant in z. □ 

3. Parabolic Schauder-type estimates 

The parabolic Schauder estimates here refer to boundedness property of the heat operator 
with respect to Holder spaces, defined in the introductory Section 1. The next auxiliary 
proposition identifies the Holder space A 2+a,1+ 2"(M x [0, T]) as a Banach space. 

Proposition 3.1. Let X be a compact manifold with boundary, C^°(X) smooth functions 
on X with compact support in the interior and &!'(X) the space of distributions. Assume 
H C &{X) is a Banach space with norm || ■ ||# such that C™(X) C H. Let D = (Di, ...,D r ) 
be a collection of differential operators on X , acting on elements of H in the distributional 



YAMABE FLOW ON MANIFOLDS WITH EDGES 13 

sense. Define 

H D := {u G H | D ± u, ..., D r u G H}, 

r 
\\u\\d '■= \\u\\h + > II-DiWIIh, U E H. 

' 4 J 

i=i 
Then (Hp, \\ ■ \\o) is a Banach space. 

Proof. Let (w n )neN C H D be a Cauchy sequence with respect to || • \\d- By completeness 
of (H, || • \\h), the sequences (u n ) n€ ^ and (DjU n ) ne ^,j = l,...,r converge to u G H and 
Vj E H,j = 1, ...,r, respectively. For any G C^(X),u G -£/£> and j = 1, ...,r we compute 

(DjU, 4>) = (u, L>*0) = lim (u n , L>*0) = lim (DjU n , <p) = (vj, <p). 

J n—KX> J n—>-oo 

Hence, DjU = Vj G H in the distributional sense. Consequently, u G H D and therefore 
(H D , || • \\ D ) is a Banach space. □ 

Since the Holder space A a '2"(M x [0, T]) and any collection of differential operators satisfies 
the conditions of Proposition 3.1, the associated spaces A 2k+Q,k+ 2'(M x [0, T]),fc G N are 
Banach spaces as well. Observe that by local parabolic regularity away from x = 0, these 
spaces coincide with classical Holder spaces in the interior, in addition to capturing the 
boundary behaviour we need for the analysis of the Yamabe flow. 

We note that by using the product rule it is easy to see that A 2k+a ' k+ 2"(M x [0, T]) is closed 
under multiplication of functions. We will also weight the spaces A 2k+a ' k+ 2'(M x [0, T]) with 
powers of t, note that u G t s A 2k+a,k+ 2"(M x [0,T]) if and only if u = t 5 u for an element 
ue A 2k+a - k+ f(Mx [0,T]). 

Theorem 3.2. Let (M m ,g) be a feasible edge space with edge data (B b ,F-F,<j)) and a singular 
neighbourhood <$/ C M of the edge B. The incomplete edge metric is given by g = go + h 
with 

9o\w = dx 2 + <f>*§ + x 2 k. 

Let Aq > / be the smallest non-zero eigenvalue of A K t,,b G B. Put 



«o := ~ h Y 7 h A - 1 > 0. 

Denote by A g the Friedrichs extension of the Laplacian on (M,g). Then for a G (0, a ) 
the heat operator e~ iAs acts as a bounded convolution operator 

e -tA g :A a,f( Mx [0,T])^A 2+Q ' 1+ f(Mx [0,T]). 

Proof. The proof we present is modeled on the classical estimation given in [LSU67]. Given 
/ G A Q, 2"(M x [0,T]), we use the asymptotics of the heat kernel to estimate each part of 
the Holder norm of e~ tAs) * f. We break this into a number of steps, performing the more 
complicated estimations first, since they introduce techniques that are used throughout the 
proof. By freezing the time variable we first estimate the Holder seminorms in space, then by 
freezing the space variables we estimate the Holder seminorms in time. Note that in each of 
these steps we will have to localize and consider arguments near each corner of the blow-up 
heat space. We finish by estimating the supremum norms at the end of the proof. 
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We will make frequent use of the fact that the edge heat kernel is stochastically complete, 
i.e. satisfies 

/ H(t,p,p)dvo\ g (p) = 1, for allp G M,t > 0, (3.1) 

J At 

which allows us to introduce Holder differences when needed. This is a simple consequence 
of uniqueness of solutions to the heat equation. More precisely, with the Friedrichs extension 
A g being a self-adjoint unbounded operator in the Hilbert space L 2 (M, vol(g)), the solutions 
u(t, •) G @{A g ) to the initial value problem 

d t u + A g u = 0, w(0) = u e 9{A g ) c L 2 (M, vol(sO), 

are unique and in fact given by u{t) = e~ tAg uo G 3t(A g ) for any t > 0. Consequently, 
it = 1 is the unique solution to the heat equation with the initial value 1 G S>(A g ), given in 
terms of the heat operator by u = 1 = e~ tAa l. This is precisely the statement of stochastic 
completeness (3.1). 

Finally, we remark that the main idea of the proof in each coordinate system is to compare 
the asymptotics of the differentiated heat kernel with the powers of defining functions that 
arise from pulling the volume form back in coordinates. The most difficult estimations that 
arise will leave an apparently singular factor of a defining function which is then appropriately 
estimated. 

Estimation of the Holder difference in space 

Consider now for any X G {A g , x^ 1 V c }, / G A Q '2"(M x [0, T]) and any p,p' G M 
Xe- tA **f(t,p)-Xe- tA **f(t,p>) 

[ XH(t-t,p,p) [f(t,p)-f(t,p)] dvol g (p)dt 

Jm+ 

[ XH(t - t, p', p) [f(t, p)-f(t, p')] dvol g Q5) dt 
Jm+ 

+ [ f [XH(t-t,p,p)-XH(t-t,p',p)] [f(t,p)-f(t,p)] dvol g (p)dt 

JO J M~ 

-ft XH(t-t,p>,p)[f(t,p)-f(t,p')]dvol g (p)dt, 

JO J M~ 

where we used (3.1) and introduced the following notation 

M + := {p G M | d M (p,p) < 2d M (p,p')}, 

M~ :={peM\ d M {p,p) > 2d M (p,p')}. 

Below we will lift the densities to the blowup space -#^. We will identify the integration 
regions M, M + and M~ with their corresponding lifts without further comment. We denote 
the four integrals above by I\,l2,l3 and I± in the order of their appearance and estimate 
these integrals separately. Note that in case of X = d t , we can employ the heat equation 

d t e- tA * * / = -A g (e- tA ° */) + /. 

Hence without loss of generality we can consider X G {A g , x~ 1 V e }. 

Estimation of the first and second integrals I±, J2. 
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We lift the heat kernel to a polyhomogeneous distribution /3*H on the parabolic blowup 
of the heat space J^%. The estimates in the interior of ^^ are classical and hence we may 
assume that /3*H is compactly supported in an open neighbourhood of the front face and 
estimate I\ and I2 near the various corners of ff. For any X G {A 9 , ;r _1 V e }, we find by 
Proposition 2.3 and (v) of Definition 1.2 

p{XH) = ( Pmd y m - 2 p$G (3.2) 

where G is a bounded polyhomogeneous distributions on ^j 1 , and by (iv) of Definition 1.2 
its th order coefficient in the right face expansion is constant in the first fibre variable. We 
emphasize that this is the worst case estimate for X = A g and that for X G {x^ 1 V e } the 
front face and temporal diagonal asymptotics of (3*(XH) is at least one order higher. 

We write down the estimates for I\. The second integral I2 is estimated along the same 
lines. 

Estimates near the lower left corner of the front face: Let us assume that the heat 
kernel H is compactly supported near the lower left corner of the front face. Its asymptotic 
behaviour is appropriately described in the following projective coordinates 

t — t x y — y ~ 

r = — — , s = -, u = , x, y, z, z, 

XZj nn nn 

JU JU 

where in these coordinates r, s, x are the defining functions of tf, If and ff respectively. 
The coordinates are valid whenever (r, s) are bounded as (t — t, x, x) approach zero. For the 
transformation rule of the volume form we compute 

(3*(dtdvo\ g (x,y,^)) = h ■ x m+2 s* dr ds du dz, 

where h is a bounded distribution on ^#^. Hence, using (3.2) we arrive for any / e 
A a, 2"(M x [0, T]) after cancellations at the estimate 

\h\<C\\f\\ a d M ((x,y,z),(x',y',z')r. 

Estimates near the lower right corner of the front face: Let us assume that the heat 
kernel H is compactly supported near the lower right corner of the front face. Its asymptotic 
behaviour is appropriately described in the following projective coordinates 

t — t x y — y _ „ 

r = -~-, s = -, u= -^r-, z,x, y, z, 

X A X X 

where in these coordinates r, s, x are the defining functions of tf, rf and ff respectively. 
The coordinates are valid whenever (r, s) are bounded as (t — t, x, x) approach zero. For the 
transformation rule of the volume form we compute 

/3*(dt dvol g (x, y, F)) = h ■ x m+l dr dx du dz, 

where h is a bounded distribution on hA\. Hence, using (3.2) we find for any / G 
A Q, f (M x [0,T]) after cancellations 



I A I < C||/IU / / x 1 dM((x,y, z),(x,y,z)) a G dr dxdudz 
J Jm+ 

The integration region lies within M + fl {x < x}, since we assume that the heat kernel 
is compactly supported near the lower right corner of the front face with s < 1. Put d := 
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dM{{x,y, z), (x',y', z')). Then for any x G M + fl {x < x} we find \x — x\ < d and x < x. 
Consequently 

x <x < x + d = x + dAf((x, y, z), (x', y', z')). 
Hence 

p px-\-d 

\h\ < c ||/|| a / / x~ 1 dM((x,y, z),(x,y,z)) a G drdxdudz 

<c'\\f\\ a / J" 1+a (|s-l| 2 + |s + l| 2 |^-z| + |u| 2 ) f Gdrdxdudz 

< c" \\f\\ a [(d M ((x iyi z),(x\y',z')) + x) a - x a ] 
= c" ll/IU rf M ((x, y, z), (*', y / , z')) f ((1 + x/rf)° - (x/rf) Q ) . 
The function w(r) = (1 + r) a — r a , r G IR + is bounded as r — y oo and hence we obtain 

\Ii\<C\\f\\ a d M ((x,y,z),(x',y',z')r. 



Estimates near the top corner of the front face: Let us assume that the heat kernel 
H is compactly supported near the top corner of the front face. Its asymptotic behaviour is 
appropriately described in the following projective coordinates 

i — ~ . x ~ x y — y ~ 

P=Vt-t, 4 = -, 4 = -, u = , y, z, z, 

P P P 

where in these coordinates p, £, £ are the defining functions of the faces ff, If and rf respec- 
tively. The coordinates are valid whenever (p, £, ^) are bounded as (t — t, x, to) approach zero. 
For the transformation rule of the volume form we compute 

P*(dt dvo\ g (x, y,z)) = h- p m+1 I s dp d£du dz. 

where h is a bounded distribution on Ki\. Hence, using (3.2) we find for any / G 
A Q, f (M x [0,T]) after cancellations 




A/+ 



Ii|<c||/|U/ / p- 1+a VU-0 2 + (£ + 2 k-21 2 +M 2 i f Gdpdidudz 



<c'\\f\\a / p- 1+a Z f dpdt=<f\\f\\ a / X f p- 2 - f+a d P dx. 

J J M+ J J M+ 

We integrate within {p > x}, since ^ < 1 near the top corner. Moreover, over M + we 
clearly have \x — x\ < d,M((%, y, z ), (x', y', z')) =: d. Hence we find 



oo 



|/i| < c" H/IU / & / p- 2 - f+a dpdx < c" ||/|| Q / x- 1+a dx 

J \x— x\<d Jx J \x—x\<d 

We consider the cases x > x and x < x separately and find 

f C" \\f\\ a d M {{x, y, z), (x', y\ z')T{{\ + x/d) a - (x/d) a ), iix<x, 
lll -\ c'" \\f\\ a d M ((x, y, z), (x', y\ z')) a {{x/d) a - (x/d - 1) Q ), if x > x. 

The functions voir) = (1 + r) a — r a and w(r) = r a — (r — l) a are both bounded as r — y oo 
and hence in both cases we deduce 

\h\<C\\f\\ a d M ((x,y,z),(x',y',z')r. 
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Estimates where the diagonal meets the front face: We proceed using projective 
coordinates (t,s,u) near the left corner, which are valid near td as well, r\ := sfr is the 
defining function of td and x the defining function of ff. For the transformation rule of the 
volume form we compute 

f3*(dt dvol g (x, y,z)) = h-7] x m+2 s^dr] ds du dz, 

where h is a bounded distribution on M.\. Hence, using (3.2) we find for any / e 
A Q, f (M x [0,T]) after cancellations 

| ii| < c||/|| Q / T]~ m ~ 1 dM{{x,y, z),(x,y,z)) a G dr/dsdudz. 
Jm+ 

We estimate the distance 



d M {{x,y,z),{x,y,z)) = y/\x - x\ 2 + (x + x) 2 \z - z\ 2 + \y - y\ 2 

= X y/\l - S\ 2 + (1 + S) 2 \Z - Z\ 2 + |«| 2 



> x y |1 — s| 2 + \z — z\ 2 + |m| 2 =: xr(s, z, z, u). 
On the other hand 



d M ((x, y, z), (x, y,z)) = x a/|1 - s| 2 + (1 + s) 2 \z - z\ 2 + 
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ur 



< 2x a/|1 — s| 2 + \z — z\ 2 + |w| 2 = 2xr(s, z, z,u). 

We transform (s, u, T) around (1, 0, z) to polar coordinates with r as the radial coordinate. 
After integration in the angular coordinates and substituting a = rj/r, we have 

|ii| < c"\\f\\ a x a f f a - m - l r- 1+a G da dr 

//•oo l"-"J 

r- l+a dr I o- m - l Gdo, 

where the integration region in r lies within {r < du {{x, y, z), (x', y', z'))/x}. a is given in 
terms of the projective coordinates (S, U, Z) near td, see (2.4) as follows 



a 



V \ 



2 / ~\ 2 /~ ^ 2 

x — x \ ( y — y \ I x(z — z) 



+ 4^4 + \-^=±\ =vWWT^. 



t-tj wt-tj \ vt-t 



Consequently the a integral in (3.3) is bounded, since G is vanishing to infinite order as 
\(S, U, Z)\ — >■ oo. As the integration region in r is within {r < d M ((x, y, z), (x 1 , y', z'))/x}, 

\h\<c"\\f\\ a x a fr- 1+a dr<C\\f\\ a dM{{x,y,z),{x',y',z')) a . 



Estimation of the fourth integral I4. 
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First we employ stochastic completeness of the heat kernel and compute 
h = [ [ XH(t - % p', p) [f(t, p) - f(t, p')] dvol g (p) dt 

= [ [f(t,p)-f(%p')] [ XH(t-t,p',p)dvo\ g (p)dt 

JO J M 

- [ [f(t,p)-f(t,p')] [ XH(t-t,p',p)dvol g (p)dt 

Jo Jm+ 

= [ [f{l P) ~ f(% p')] [ -XH(t - I p\ p) dvol g (p) dt. 

Jo J M+ 

We lift the heat kernel to a polyhomogeneous distribution j3*H on the parabolic blowup 
of the heat space ^^, and as before assume that (3*H is compactly supported in an open 
neighbourhood of the front face and estimate I4 near the various corners of ff. Estimation 
near the lower left corner of the front face follows along the lines of I\ estimates. The 
estimates near the top corner are parallel to those near the right corner. Therefore, we only 
explicate the estimates near the right corner and the diagonal. 

Estimates near the lower right corner of the front face: 

The asymptotic behaviour of (3*H is appropriately described in the following projective 
coordinates 

t — t x' y' ~y 1 ~ ~ 

t = -=«-, s = —, u= - , z ,x, y, z, 

where in these coordinates r, s, x are the defining functions of tf, rf and ff respectively. 
The coordinates are valid whenever (r, s) are bounded as (t — t, x' , x) approach zero. For the 
transformation rule of the volume form we compute 

(3*(dtdvol g (x, y, z)) = h ■ x m+1 dr dx du dz, 

where h is a bounded distribution on J\A\. We separate X into two components X = 
X' + X", where all components which contain derivatives in y are comprised into X", so that 
we have the following worst case estimates 

(3*{X'H) =x- m - 2 s ao T°°G, 

(3*(X"H) = J-^Vt 00 ^. 

We separate J 4 = I4+/4 correspondingly. The estimation of I' A is easy, since by Proposition 
2.3 the kernel (3*(X'H) is at least of leading order ao at the right face. We find for any 
/ G A Q, f (M x [0,T]) after cancellations 

l^l ^ c\\f\\ a d M ((x,y, z),(x',y\ z')) a / x _1 s a ° G dr dx du dz 

Jm+ 

/■oo 



< c \\f\\ a d M ((x, y, z), (x', y', z')) a x^ J x~ x ~^dx 

<d'\\f\\ a d M {{x,y,z)\x\y\z')Y. 

We exemplify estimation of I'[ for X" = d V:j o X'" where yj eRisa local coordinate on B, 
and perform integration by parts in the corresponding Uj G M. For fixed (x, z) the integration 
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region is Uj G [— R(x, z), R(x, z)], where \uj\ = R(x, z) corresponds to (x, y, z) G dM + . Write 

, ..,Uj, ..,' 



u := (ui, ..,Uj, ..,Ub) G M} 1 . Integration by parts gives 



I'l = I Sf [ x- l d u {f3*X'"H) x m+1 h(x, y-xu,z) dr dxdu dz 
J Jm+ 

= I Sf fx- 1 f3*X'"Hx m+1 h(x,y-xu,T) drdxdudz 

5f J x- 1 p*X' n HW m + x d VL h(x, y-xu,z) drdxdu dz 

Jm+ 

= A + B, 

where we simplified the expression with 

Sf ■= fit - rx 2 , x, y, z) - fit - tx 2 , x', y', z'). 

Since h is smooth in the edge variable, d u h(x,y — xu, z) = xd y h(x,y — xu, z) and we deduce 
after cancellations 

\B\<c\\f\\ a d M ((x,y,z),(x',y',z')r. 

For the first integral we note that for \v,j\ = R(x, z) we have (x, y, ~z) G dM + so that 
d M ((x, y, z), (x, y, J)) = 2d M ((x, V, z), (a/, y', z')), 
dAi((x', V, z'), (x, y, z)) > d M ((x, y, z), (x, y\ z')), 



where the second relation follows from the first one using the triangle inequality. Using 
the second relation we find 



\A\<c\\f\\ a d M ((x,y,z),(x',y',z')r Jx-'G 



dr dx du dz 

\uj\=R(x,z) 



^ C II/IU / x 1 d>M{{x' -,y' i z'),(x,y,z)) a G dr dxdudz 



Proceeding exactly as for the estimate of I\ at the right face, we obtain 

\A\<c\\f\\ a d M ((x,y,z),(x',y',z')r. 

Estimates where the diagonal meets the front face: 

The asymptotics of the heat kernel near the diagonal meeting the front face are conveniently 
described using the following projective coordinates 

,f = tzl, S = %=R, U = 4=4, Z = ^#, A y>, J. (3.4) 



X' 



't-t y/t-t' y/t-t 



In these coordinates tf is the face in the limit \(S, U, Z)\ — > oo, ff and td are defined by 
x', rj, respectively. For the transformation rule of the volume form we compute 

(3*{dtdvol g {x,y,T)) = h ■ (x') m+ V +1 (l - r]S) f d7]dSdUdZ, 

where h is a bounded distribution on hA\. Hence, using (3.2) we find for any / G 
A Q 'f (M x [0,T]) after cancellations 

|/ 4 | < c\\f\\ a d M {{x,y,z),{x',y',z')) a f f rj^hGdridSdU dZ. 

J Jm+ 
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We see that the crude estimate (3.2) leads to an a priori singular behaviour in rj. In order 
to overcome this difficulty note first the following transformation rules for the edge vector 
fields V e 

(3*(x'd x/ ) = -2 v d v + -d s + Zd z + x'd xl , f3*(x'd y ,) = -du + x'dy, 0*(d z ,) = -d z + d z ,. 

7] rj rj 

We find that the singular behaviour in rj comes from differentiation in (S,U,Z). Hence 
we might as well consider the case j3*(XH) = (x'r])~ 2 dgP*H. Estimates of the other cases 
follow along the same lines. For fixed (rj, U, Z) the integration region is {\S\ < R(i], U, Z)}, 
where \S\ = R(r), U, Z) = R corresponds to (x, y, z) G dM + . Integration by parts gives 



Sf / (x'ri)- 2 d 2 s (P*H)h(x') m+2 V m+1 (i ~ vS) f dr]dSdUdZ 

Jm+ 

Sf f (x'ri)- 2 d s ((3*H)h(x') m+2 ri m+1 (l - r]S) f 
Jm+ 



M+ 

dr] dU dZ 



\S\=R(ri,U,Z) 

Sf f (x'ri)- 2 d s (f3*H)d s (h(l - r]Sy)(x') m+2 r] m+1 dr]dU dZ 

JM+ 

— V -t- T" 

where we have simplified the expression with 

Sf := fit - (x'r)) 2 , x, y, z) - f(t - (x' V ) 2 , x', y', z'). 

The asymptotic behaviour of the heat kernel implies ds(P*H) = (x'r])~ m( ^, where & is 
bounded and vanishes to infinite order as 1(5, U, Z)\ —¥ 00. Using the fact that h is a smooth 
function of x = x'(l—r]S) we find after cancellations, that I'[ does not have singular behaviour 
in rj and estimate 

\I';\<C\\f\\ a d M ((x,y,z),(x',y' } z')r. 

For the estimation of I' A , we note that for |,S| = R(i], U, Z) we have (x, y, F) G dM + so that 
d M ((x, y, z), (x, y,z)) = 2d M ((x, y, z), (x', y' , z')), 
d M ((x', y\ z'), (x, y, z)) > d M ((x, y, z), (x', y', z')), 
d M ((x', y', z'), (x, y, z)) < 3d M ((x, y, z), (2', y\ z')), 



where the second and third relations follow from the first one using the triangle inequality. 
We now transform to coordinates which are valid near the left corner of the front face and 
also near the diagonal 

2 _ t-t _x_ _ y' -y , , , ~ 
V / ,191 s ,1 u , j x , y , z , z, 

{x') z x' x' 

where in these coordinates (r], s, x') are the defining functions of tf, If and ff respectively. 
The coordinates are valid whenever (rj, s) are bounded as (t — t, x', x) approach zero. The 
volume dU dZ transforms as follows 

dUdZ = ri~ m+1 dudz. 
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Hence we find 

\I' 4 \<c\\f\\ a d M ((x,y,z),(x',y',z')) a \ rf m <S d v dudl 

J \S\=R 

drj du dz 

\S\=R 

Now we proceed exactly as before in the estimation of I\. We pass from (w, z) to polar 
coordinates around (0, z') with the radial function 



<c'||/|U V- m &d M ((x',y',z'),(x,z,z)y 



r(u,z',z) = a/I-z' — z\ 2 + \u\ 2 . 
Observe the following relation 



c < d M ((x',y',z'),(x,z,I)) < c/ y /\S\* + \U\i+]Zf 
x'r{u,z',z) ~ ^\U\ 2 + \Z\ 2 

Consequently substituting a = rj/r and using r < d((x,y,z),(x,y,z))/x', which follows 
easily from the first inequality above, gives 



\I'A<c"\\f\L(xT I I r - 1+a a- m 



^\S\ 2 + \U\ 2 + \Z\ 2 



o J VW< 

<c'"\\f\\ a d((x,y,z),(x,y,I)r. 



dr da 
\s\=R 



Estimation of the third integral I3. 

As in the previous estimates we lift the heat kernel to a polyhomogeneous distribution /3*H 
on the parabolic blowup of the heat space ^#^, and assume that /3*H is compactly supported 
in an open neighbourhood of the front face. Using the standard coordinates and the mean 
value theorem we obtain 

t 



I 3 = \x- x'\ / / d(.XH{t - t, £, y, z, x, y, z) [fit, x,y,z) - f(t, x, y, z)] dvol g dt 

JO J M- 

+ \y-y'\ / d^XHit - t, x', 7, z, x, y, z) [f(t, x,y,z) - fit, x, y, z)] dvol g dt 
Jo Jm- 

+ \z- z'\ / d c XH{t - t, x', y', C, x, y, z) [f(t, x,y,z)- f(t, x, y, z)] dvol g dt, 

Jo J M- 

for some intermediate values (£,7, C)- Consider any (x",y",z") G M with 

d M (x, y, z, x", y", z") < d M (x, y, z, x' , y', z'). 
Then for any (x, y,T) G M~ we compute 

2d M (x, y, z, x', y', z') < d M (x, y, z, x, y, z) 

< d M {x,y,z,x ,y ,z ) + d M {x ,y ,z ,x,y,z) 

< d M (x, y, z, x', y , z') + d M (x", y", z", x,y,z). 
We deduce 
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d M (x,y,z,x',y',z') <d M (x",y",z",x,y,z), 
d M (x, y, z, x,y,z)< 2d M (x", y", z", x,y,z). 

Using these inequalities we may estimate ^3 for any / e A a '2"(M x [0, T]) as follows 

\h\ < 2 \\f\U\x - x'\ \ \ \d(X.H{t - t,£,y, z,x,y, z)\d M (£,y, z,x,y, z) a dvol g dt 
Jo Jm- 

+ 2\\f\\ a \y-y'\ / / \d 1 XH(t-t,x',j,z,x,y,z)\d M (x',"f, z,x,y,z) a dvol g dt 
Jo Jm- 

+ 2\\f\\ a \z-z'\ [ [ \d c XH(t-lx',y',C,x,y,I)\d M (x',y',C,x,y,I) a dvol g dt 
Jo Jm- 

Estimates near the lower left corner of the front face: 

Let us assume that the integrand in each of the components of ^3 is compactly supported 
near the lower left corner of the front face. Their asymptotic behaviour is appropriately 
described in the following projective coordinates 



t — t x y — y 



\h) T = —£Ti s = 7' u = —^i * = S> V, z > z i 



( J 3) r = 7-A2> s = — ' u = — — ' e = x '> y> z > *> ( 3 - 6 ) 

(x ) l X X 

(this t-t x y'-y , 

[x'Y x' x' 

where in these coordinates r, s, 9 are the defining functions of tf, rf and ff respectively. 
The coordinates are valid whenever (r, s) are bounded as (t — t, 9, x) approach zero. For the 
transformation rule of the volume form we compute 

/3* (dt dvolg (x, y,z)) — h- 9 m+2 dr ds du dz, 

where h is a bounded distribution on ^#^. Hence, using (3.2), and regarding (x, y) as 
functions of (s,u), we arrive for any / e A a '2"(M x [0, T]) after cancellations at the estimate 
(note, differentiation in ( does not lower the front face asymptotics) 

|l£|<c||/|| a |x-x / |r 1+a / / (d M (^,y,z,x,y,T)/O a GdTdsdudz, 

J Jm- 

I/3I <c\\f\\ a \y-y'\(x')- 1+a J J (d M (x', 1 ,z,x,y,T)/x') a GdrdsdudI, 

l/3'l <c\\f\\ a \z-z'\(x') a f f (d M (x',y',(,x,y,7)/x') a GdTdsdudz. 

J Jm- 
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Note that in each case we can estimate the integrand using the appropriate projective 
coordinates 



(I3) d M (£, V, z, x, y, z)/£ < 2a/|1 - s\ 2 + \z - z\ + \u\ 2 

(Jg) d M (x', 7, z, x, y, z)/x' < 2a/|1 - s\ 2 + \z-z\ + \u\ 2 \ =: 2r(s, u,z). 

(if) d M (x',y',(,x,y,z)/x' < 2^1 - s\ 2 + \(-z\ + W 2 

By (3.5) the integration region M~ lies within {r > cIm(x, y, z, x', y', z')/9}. If u is bounded, 
then r is bounded near the left corner, so that 1 < cr~ x and hence we may estimate constants 
against as many negative powers of r as we like. Similarly, if \u\ tends to infinity, then r/\u\ 
is bounded near the left corner, and we may introduce as many negative powers of r/\u\ as 
we like. The additional powers of \u\ are irrelevant, since G vanishes to infinite order as 
|w| —7- 00. Hence we may estimate 

%\ <d\\f\\ a \x- x'\ r 1+a (d M (x, y, z, x', y', z')/0~ 1+a , 

\^\<c'\\f\\ a \ y -y'\(xr 1+a (dM(x,y,z,x',y',z')/xr 1+a , 

\I^'\ < c' \\f\\ a \z - z'\(xT(dM(x,y,z,x',y',z')/xT 1+a - 
Altogether we finally obtain 

\h\ < c" ll/IU (\x - x'\ + \y - y'\ +x'\z-z'\) d M (x, y, z, x', y', z')- 1+a 

<c'"\\f\\ a d M (x,y,z,x',y',zT- 

Estimates near the lower right corner of the front face: 

Let us assume that the integrand in each of the components of I3 is compactly supported 
near the lower right corner of the front face. Their asymptotic behaviour is appropriately 
described in the following projective coordinates 

(i 3 ) t = —-, s = -, u = —=—, 9 = f , x, y, z, z, 

JU JLj ju 

/r/A t-t x' j-y , „ 

(i 3 ) r = — -, s = — , u = —=—, 9 = x, x, y, z, z, 

JU iaj JU 

/w/a t-t x' y'-y , „ 

(i 3 ) t= -— , s = — , u = -^—, 9 = x , x, y, C, z, 

X z X X 

where in these coordinates r, s, x are the defining functions of tf, If and ff respectively. 
The coordinates are valid whenever (r, s) are bounded as (t — t, 9, x) approach zero. For the 
transformation rule of the volume form we compute 

/3*(dt dvol g (x, y,'z)) — h- x m+l dr dxdu dl, 

where h is a bounded distribution on ^ 2 . Hence, using (3.2) we arrive for any / G 
A a '2"(M x [0,T]) after cancellations at the estimate (note, differentiation in £ does not lower 
the front face asymptotics) 

\I' 3 \ < c\\f\\ a \x - x'\ j I x- 2+a s~ 1+ao (d M (t,y,z,x,y,z)/x) a GdTdxdudz, 
l-JJs'l < c||/|| tt |j/-|/| / / x' 2+a (d M (x',-f,z,x,y,z)/x) a Gdrdxdudz, 
\l'a\<c\\f\\ a \z-z'\ f f x- 1+a s ao (d M (x',y',(,x,y,Z)/x) a GdTdxdudI. 
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The additional factors s~ 1+a ° and s a ° require explanation. Proposition 2.3 and condition 
(iv) of Definition 1.2 imply together that near the right corner as s — > the leading s° 
coefficient in the expansion of (3*(XH) is fibrewise constant and hence vanishes under differ- 
entiation in ( and £. By construction, the next term in the expansion of fi*{XH) is of the 
form a ■ s a ° , or higher. 

We begin with the estimation of Jg. Put d := dM(x,y,z,x',y',z'). We separate the 
integration region into M{~ = M~ fl {x < d} and M^ = M~ D {x > d}. Then using 
\y — y'\ < d and (3.5) we find 



|/?l < c' 



+ c' 



<c" 



+ c" 



a\y-y 



\y-y 




Mf 




x +a (dM(x', 7, z, x, y, z)/x) a G dr dx du dz 



x " +a ((iM(a ;/ ,7, z^x^y^/x) 01 G dr dxdudz 



M~ 




x +a (dM{x' ,7, z, x,y, z)/x) +a Gdrdxdudz 



Aff 




x 2+q (o?m(^ / ,7, z,x, y,z)/x) a Gdrdxdudz. 



m: 



Note that using the appropriate projective coordinates we find 



(4) d M(£,y,z } x,y,z)/x < 2^(1 — s| 2 



(/g) d M (a/, 7, «, », ?/, F)/x < 2a/J1 
(Jg") ^ / (x', 2 /',C,5?,^i)/^<2 v ^T 3 



+ |^r — 5] + 



si 2 + \z — 



IC 




=: 2r(s, u, z). 



(3.7) 



If u is bounded, then r is bounded near the right corner. If \u\ tends to infinity, then r/\u\ 
is bounded near the right corner, and the additional powers of \u\ get absorbed by G which 
vanishes to infinite order as \u\ — > oo. Hence we can estimate J rG du against a constant and 
find 



|/ 3 '| < c'" 



X 



-1+a 



dx 



+ c" 



d 



X 



-2+a 



dx 



<c 



d° 



In order to estimate J 3 ", note that sx = x' and hence 



l^ / l<t/||/|L(x , ) ao U-z , | 



X 



-l-^o+a d ~ 



<c" 



< c" \\f\\ a (x + x') ao \z - z'\ a < C \\f\\ a d a . 



\z-z'\(x') a 



It remains to estimate 7 3 . Here we employ the more refined statement of Proposition 2.3 on 
the front face behavior of coefficients. The proposition asserts that H = H' + H", where /3*H' 
is of leading order (— m) at the front face and has only terms of the form s 1 in its expansion 
at the right face, where 7 are the indicial roots of A g . f3*H" is of leading order (— m + 1) at 
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ff and its expansion at the right face is comprised of terms s 7+/ , I > 1. Consequently 

1-^3 1 < c \\f\\a \x — x'\ I I d^(5*{XH'){dM{i ) y ) z,x,y,z)/x) a dr dxdudz 
J Jm- 

+ c \\f\\a\x — x'\ / / d^f3*(XH")(dM(^,y, z,x,y,z)/x) a dr dxdudz 
— c' \\f\\ a \x — x'\ j I x^ 2+a s a °{d M {i,y,z,x,y,z^)/x) a G dr dxdud'z 

+ c' \\f\\ a \x-x'\ ff x- 1+a s- 1+ao (d M (^y,z,x,y,I)/x) a GdTdxdudI. 

The first integral is estimated precisely as 1'^ . For the second integral we find (note sx = £) 
using (3.5) 

c ' \\f\\a \x — x'\ / / x~ 1+a s~ 1+a °(dM{£,y,z,x,y,T)/x) a Gdr dxdudz 

|x ~ x| j / / x- ao+a {d M {£,y,z,x,y,z)/x) a GdT dxdudz 

<c"\\f\\ a \x-x'\ ao I I x- ao+a (d M (t,y,z,x,y,z)/x) 1 - ao+a GdTdxdudz' 

J J M~ 

<c"||/|U|a;-a;r o / x~ ao+a dx < C \\f\\ a d M (x,y, z, x',y', z') a - 
Jo 

Estimates near the top corner of the front face: Let us assume that the heat kernel 

H is compactly supported near the top corner of the front face. Its asymptotic behaviour is 

appropriately described in the following projective coordinates 



(h) P = V t - t, v = -, v = -, u = , 9 = f , y, z, z, 

P P P 



j 



[I'D p = y/t-t, v = — , v = -, u = , 9 = x', 7, z, z, 

P P P 



x' _ x y' — y 



{I'i') P = Vt-t, v = — , v = -, u = , 9 = x', y', C, z, 

P P P 

where in these coordinates p, u, v are the defining functions of the faces ff, If and rf respec- 
tively. The coordinates are valid whenever (p, v, v) are bounded as (t — t, x, £, x', x) approach 
zero. For the transformation rule of the volume form we compute 

f3* (dt dvolg (x,y,z)) = h ■ p m+1 W dp dv du dz. 

where h is a bounded distribution on M. 2 h . Hence, using (3.2) we find for any / e 
A Q, f (M x [0,T]) after cancellations 



\I',\<c\\f\\ a \x-x'\ / / p- 2+a v- 1+ao (d M (Z,y,z,x,y,z^/xrGdpdvdudz, 
J Jm- 

\I'z\<c\\f\\ a \y-y'\ / / p~ 2+a (d M (x',~f,z,x,y,z)/x) a Gdpdvdudz, 

J J M~ 

14" I < c\\f\\ a \z - z'\ II p- 1+a u ao (d M (x',y',C,x,y,I)/x) a GdpdududI. 
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The additional factors i/~ 1+a ° and v a ° arise as before in the discussion near the right corner. 
The integrals are estimated in exactly the same manner as near the right corner. 

Estimates near the diagonal meets the front face: We proceed using projective 
coordinates (3.6) near the left corner, which are valid near td as well, rj := y/r is the defining 
function of td and x the defining function of ff. For the transformation rule of the volume 
form we compute 

/3*(dt dvol g (x, y,z)) = h-T] 8 m+2 s^drj ds du dz, 

where h is a bounded distribution on M.\. Hence, using (3.2) we find for any / e 
A Q, 2'(M x [0, T]) after cancellations (note, differentiation in ( does not lower the front face 
asymptotics) 



141 


<c||/||a 


\x 


1 
— X 


ir-j 


141 


<c||/|| Q 


\y 


-y'\ 


(x'y 1+a 


41 


<c||/|U 


\z 


-A 


(x'r [ 




V m (dM(£,,y, z,x,y,~z)/^) a G drjds dudz, 



M~ 




i] m {dM{x' ,7, z,x,y,z)/x') a G drjds dudz, 



M- 



V m (dM(x',y',C,x } y,z)/x') a G drjds dudz. 

J J M- 

We estimate the distance term using appropriate projective coordinates 

(I' 3 ) d M (Z,y,z,x,y,z)/£ < 2^1 - s| 2 + \z-z\ + |u| 2 

(40 d M {x', 7, z, x, y, z)/x' < 2 v / |l-s| 2 + |,2-z| + H 2 ) =: 2r(s, u,z). 



{I'») d M (x',y'X,x,y^/x' < 2^1 - s| 2 + |C - *| + \u\\ 

We transform (s, u, ~z) around (1,0, z) to polar coordinates with r as the radial variable. 
After integration in the angular variables and substituting a = rj/r we have 

141 <c||/|U|x-x'|r 1+a / / a- m ~ 2 r- 2+a G da dr, 

J J M- 
i\ f^>\-l+a / / „-m-2—2+a 



141 <c\\f\\ a \y-y'\{x')- l+a / / a- m - 2 r- 2+a Gdadr, 

I4'| <c\\f\\ a \z-z'\(x') a J J a- m - x r- 1+a Gdadr. 

where by (3.5) the integration region M~ lies within {r > d M (x,y,z,x',y',z')/9}. a is 
given in terms of the projective coordinates (S, U, Z) near td, see (2.4) as follows 



a' 1 = vW+MWP. 

Consequently the a integral in (3.3) is bounded, since G is vanishing to infinite order as 
| (S, U,Z)\ — > oo. We find, with the integration region in r within {r > d)tf(x, y, z, x', y', z')/9} 

141 < C 'II/IU \x-x'\d M (x,y, z,x',y',z')~ 1+a 
\I'3\<<S\\f\\a\y-y , \d M (x,y,z,x',y , ,z , )- 1+a , 

\^\<c'\\f\\ a {x')\z-z'\d M {x,y,z,x\y\z')- l+a - 
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Altogether we finally obtain 

\h\<c"\\f\\ a (\x-x , \ + \y-y , \+x'\z-z , \)dM(x,y,z,x',y , ,zT 1+a 
<c'"\\f\\ a d M (x,y,z,x',y',zT. 

Estimation of the Holder difference in space without differentiation. 

In view of the Holder norm for A 2+a ' 1+ 2"(M x [0, T]) we also need to estimate the Holder 
difference in space without differentiating the heat kernel. In other words we need to establish 
mapping properties of the heat operator as a convolution operator between A Q '2"(M x [0, T]) 
and itself. Here we use mean value theorem, (3.1) and estimating in a manner similar to I3 
above 

e- 4A * * f(t, x, y, z) - e~ tA ° * f(t, x', y', z') 

= \x-x'\ / d^H(t - t, f , y, z, x, y, z) [f(t, x,y,z) - f(t, x, y, z)] dvol g dt 
Jo Jm 

+ \y-y'\ / d~/ H (t - t, x', 7, z, x, y, z) [fit, x,y,z) - f(t, x, y, z)] dvol g dt 
Jo Jm 
ft 



+ \z- z'\ / d c H(t - t, x', y', C, x, y, z) [f(t, x,y,z) - f(t, x, y, z)] dvol g dt. 
Jo Jm 

The estimates required here are much simpler than those for 1% as we apply two fewer 

derivatives to the heat kernel. For brevity we omit the straightforward computations which 

consist only of lifting the integrand to ^£ and checking the powers of defining functions at 

the various boundary faces. Using \psPtdG\ < Cy/i we find 

||e- 4A * * f(t, x, y, z) - e-^ * f(t, x', y', z') \\ x < CV~t\\f\\ a d M ((x, y, z), (x\ y', z')) a . (3.8) 

Estimation of the Holder difference in time 

Fix t > t' without loss of generality. Assume first (2f — t) > 0. Then we consider for any 
X G {A 9 , x^ 1 V c } the difference 

Xe- tA **f(t,p)-Xe- tA **f(t',p) 

r XH(t-t,p,p) [f(t,p)-f(t,p)} dvo\ g (p)dt 



2t'-t JM 

t' 

' XH(t' -t,p,p) [f(t,p) - f(t,p)] dvo\ g (p) dt 



ll'-t JM 
2t'-t 



+ / / [XH(t-t,p,p)-XH(t'-lp,p)] [f(t,p)-f(t,p)] dvol g (p) dt 

Jo Jm 

=: Ji + J 2 + J 3 , 

where again we have used stochastic completeness of the heat kernel. Note that in case of 
X = d t , we can employ the heat equation 

d t e- tA * * / = -A 9 (e- tA * */) + /. 

Hence without loss of generality we can consider X G {A 3 , x^ 1 V e }. 
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We lift the heat kernel to a polyhomogeneous distribution /3*H on the parabolic blowup 
of the heat space J^^. The estimates in the interior of ^^ are classical and hence we may 
assume that /3*H is compactly supported in an open neighbourhood of the front face and 
estimate J\ and J 2 near the various corners of ff. By Proposition 2.3 and (v) of Definition 
1.2 we find 

0*(XH) = ( PsPtd )- m - 2 p~G (3.9) 

where G is a bounded polyhomogeneous distributions on M^. We emphasize that this is 
the worst case estimate for X = A g and that for X G x~ l V e the front face and temporal 
diagonal asymptotics of (3*(XH) is at least one order higher. 

We write down the estimates for J\ and J 3 . The second integral J 2 is estimated along the 
lines of J\. For J 3 we employ the mean value theorem and find for some 9 G (£', t) 

p2t'-t r- 

J 3 = / / [XH(t-t,p,p) - XH(t' -t,p,p)][f(t,p) - f(t,p)] dvol g (p)dt 

J0 J M 

= \t-t'\ / d e XH(0 -t,p,p) [f(t,p) - f(t,p)] dvol g (p) dt. 

JO J M 

Estimates near the lower left corner of the front face: 

Let us assume that the heat kernel H is compactly supported near the lower left corner of 
the front face. Its asymptotic behaviour is appropriately described in the following projective 
coordinates 

t — t x y — y „ 

t - , s , u , x, y, z, 2, 

x A x x 

where in these coordinates r, s, x are the defining functions of tf, If and ff respectively. 
The coordinates are valid whenever (r, s) are bounded as (t — t, x, x) approach zero. For the 
transformation rule of the volume form we compute 

(3* (dt dvo\ g (x , y, F)) = h ■ x m+2 s* dr ds du dz, 

where h is a bounded distribution on ^%. Hence, using (3.9) we arrive for any / e 
A a '2"(M x [0, T]) after cancellations at the estimate (note that near the left corner x 2 > (t—t)) 



\Ji\ < C||/||q / x a (d M ((x,y,z),(x,y,z))/x) a GdTdsdudz 

<C"||/|| a f x- 2+a ^\l - s\ 2 + \z- I\ 2 + \u\ 2 a Gdtdsdudz 

<C"\\f\\ a f {t-t)-^dt = C'"\\f\\ a {t-t'f>. 
Jit'-t 

We estimate J 3 using the projective coordinates as above where t is replaced by 6. Using 
that for t G [0, It' — t] we have \9 — t\ > \t — t'\, since we have assumed (2t' — t) > 0. Hence 
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we find 

\M <C\\f\\ a (t-t') ( x- 2+a (d M ((x,y, z), (x,y,z))/x) a Gdr ds dudz 

/■2t'-t 

< c" ii/iu (t - o / (* - *)~ 2+f dt < c"> ii/iu (t - o f • 



Estimates near the lower right corner of the front face: 

Let us assume that the heat kernel H is compactly supported near the lower right corner of 
the front face. Its asymptotic behaviour is appropriately described in the following projective 
coordinates 

t — t x y — y ___ 

t = -=5-, s = -, u= —=—, z,x, y, z, 

where in these coordinates r, s, x are the defining functions of tf, rf and ff respectively. 
The coordinates are valid whenever (r, s) are bounded as (t — t, x, x) approach zero. For the 
transformation rule of the volume form we compute 

/3*(dtdvo\ g (x, y, z)) = h ■ x m+1 dr dx du dz, 

where h is a bounded distribution on hA\. We employ the more refined statement of 
Proposition 2.3 on the front face behavior of coefficients. The proposition asserts that H = 
H' + H", where f3*H' is of leading order (— m) at the front face and has only terms of the form 
s 7 in its expansion at the right face, where 7 are the indicial roots of A g . fi*H" is of leading 
order (— m + 1) at ff and its expansion at the right face is comprised of terms s 1+l ,l > 1. 
Consequently 

\P*{XH')\ = x- m - 2 s ao G', \(3*(XH")\ = s- m_1 G". 

Hence, we find for any / e A a, a'(M x [0,T]) after cancellations 

\Ji\ <C||/|| Q f (x- 1+a s ao G' + x a G") (d M ((x,y,z),(x,y,I))/x) a dTdxdudI 

< a ii/ii. / r - s » G < + x- G ») VW^TJ^WTW S***. 

Using the fact that G vanishes to infinite order as \u\ — > oc and x > x near rf, we find 

//»oo 
(x -1 s ao G" + G") ^\l- s\ 2 + \z-z\ 2 + \u\ 2 a dxdu dz < x aa I x~ x ~^dx + c<c'. 

Consequently, using x 2 > (t — t) we find 

|^l| < C" H/lla J (x- 3+a S ao G' + X- 2+a G") y/\l-s\ 2 + \ Z -z\ 2 + \u\ 2 "dtdxdudz 

<C"\\f\\ a f {t-t)- 1+ ^dt<C'"\\f\\ a {t-t')^. 

J2t'-t 
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We estimate J3 using the projective coordinates as above where t is replaced by 9. Using 
\8 — t\ > It — t'|, we find as before 



\M <C\\f\\ a (t-t') / (x' 3+a s a °G' + x' 2+a G")(d M ((x,y,z),(x,y 7 T))/x) a drdxdudz 
<C'\\f\\ a {t- t') f (x- 5+a s ao G' + x~ 4+a G") v/|l-s| 2 + |z-z| 2 + H 2 a dtdxdu dz 

< C" ll/IU (t - / (9 - t)- 2+ Ut < C" H/lla (t - O f • 

JO 

Estimates near the top corner of the front face: 

Let us assume that the heat kernel H is compactly supported near the top corner of the 
front face. Its asymptotic behaviour is appropriately described in the following projective 
coordinates 

1 — ~ x ~ x y — y ~ 

p=Vt-t, 4 = -, 4 = -, u = , y, z, z, 

P P P 

where in these coordinates p, £, £ are the defining functions of the faces ff, If and rf respec- 
tively. The coordinates are valid whenever (p, £, ^) are bounded as (t — t, x, u) approach zero. 
For the transformation rule of the volume form we compute 

/3* (dtdvol g (x, y,z)) = h- p m+1 I s dp d£du dz. 

where h is a bounded distribution on M\. Hence, using (3.9) we find for any / e 
A Q 'f (M x [0,T]) after cancellations 

\h\ < C\\f\\ a [ p- 1+a (d M ((x,y,z),(x,y,I))/p) a Gdpd^dudI 



<C'\\f\\ a [ 2{ ' ' p- l+a dp = C'"\\f\\ a {t-t>f>. 
Jo 

We estimate J3 using the projective coordinates as above where t is replaced by 9. Using 
\9 — t\ > \t — t'|, we find as before 

\J 3 \<C\\f\\ a (t-t') [ p- 3+a (d M ((x,y,z),(x,y,I))/p) a Gdpd^dudz 



< C ii/iu (t - / p~ 3+ad P = C" ll/IU (* - *>■ 

Jy/(t-V) 

Estimates where the diagonal meets the front face: The asymptotics of the heat 
kernel near the diagonal meeting the front face is conveniently described using the following 
projective coordinates 

2 t-t ix-x) y-y x{z-z) 

v = -^-> s = -, — s= » u = -, — s= > z = -, — ^ , x > y> z - 



x y/t-t \/t-t y/t-t 

In these coordinates tf is the face in the limit \(S, U, Z)\ —$■ 00, ff and td are defined by 
x, f], respectively For the transformation rule of the volume form we compute 

0*(dtdvol g {x,y,z)) = h ■ x m+2 7] m+1 {l - r]S) f dr]dSdUdZ, 
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where h is a bounded distribution on Ai\. Hence, using (3.9) we find for any / G 
A Q, f (M x [0,T]) after cancellations 



|Ji| <C||/|| a / V - 1+a (d M ((x,y,z),(x,y,I))/ V ) a Gd V dSdUdZ 
<C"||/|| Q a; Q f V~ 1+a V\S\ 2 + W\ 2 + \Z\ 2 Gdr]dS dU dZ 

<C"\\f\\ a x a / V- 1 + a dr, = C"'\\f\\ a (t-l!)S. 







We estimate J3 using the projective coordinates as above where t is replaced by 9. Using 
\9 — t\ > \t — t'\, we find as before 



\J 3 \ <C\\f\\ a (t-t')x- 2 j r ] - 3+a (d M ((x,y,z),(x,y,I))/7 1 rGdr ] dSdUdZ 
<C'\\f\\ a (t-t')x- 2+a f 7]- 3+a ^\S\ 2 + \U\ 2 + \Z\ 2 Gdr] dS dU dZ 

POD 

< c" ii/iu (t - t')x- 2+a / v- 3+a d v = C" ii/iu (t - t') f . 



This concludes the estimation near the diagonal assuming 2t' — t > 0. 

Finally, if (2t' — t) < then we consider for any X G {A 3 , x~ x V c } the difference 

Xe- tA **f(t,p)-Xe- tA **f(t',p) 
= [ [ XH(t-t,p,p) [f(t,p)-f(t,p)] dvo\ g (p)dt 

JO J M 

- I I XH(t'-t,p,p) [f(t,p)-f(t,p)] dvo\ g (p)dt 
Jo Jm 

=:J 1 + J 2 , 

where we have used stochastic completeness of the heat kernel as before. The integrals 
J[, J 2 are estimated exactly as before the respective integrals Ji, J2, using that for (2t' — t) < 
we have t,t' <2\t-t'\. 

Estimation of the Holder difference in time without differentiation. 

As before in the parallel case of estimating the Holder difference in space without differ- 
entiation, we also need to estimate the Holder difference in time without differentiating the 
heat kernel. In other words we need to establish mapping properties of the heat operator as 
a convolution operator between A Q, 2"(M x [0, T]) and A a '2"(M x [0,T]). Here we use mean 
value theorem, (3.1) and estimate as in J 3 above (£' < t) 

e- tA * * f{t, x, y, z) - e-*' A « * f(t', x, y, z) 

= \t - t'\ / / d e H(9 - t, x, y, z, x, y, z) [f(t, x,y,z) - f(9, x, y, z)] dvol g dt. 
Jo Jm 
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The estimates required here are much simpler than those for J 3 as we estimate with two 
fewer derivatives applied to the heat kernel. For brevity we omit the straightforward compu- 
tations which consist only of lifting the integrand to ^ 2 and checking the powers of defining 
functions at the various boundary faces. Using \t — t'\ < \fi\t — t'\* we find 

||e- tA * * f(t, x, y, z) - e- 4 ' A * * /(£', x, y, z)\\ a/2 < Cy/i\\f\\ a \t - i'|l (3.10) 



Estimation of the Supremum Norm 

For any X G {A 9 , x~ x V e } stochastic completeness of the heat kernel yields 

|*e- iA **/| < H/IU / / XH(t-lp,p)d M (p,p) a dvo\ g (p)dt. 
Jo Jm 

As before, the case X = d t can be reduced to X G {A 5 , x~ x V e } using the heat equation. 
The integral above can be estimated against a constant in a straightforward ways using 
Theorem 2.2 and Proposition 2.3, in the projective coordinates near the various corners of 
the front face. The estimates near the diagonal are not so straightforward but follow the 
arguments for estimation of I\ near the diagonal along the same lines. 

In the case where X = id, we cannot use stochastic completeness of the heat kernel, but 
the estimates are straightforward. 

The proof of the Theorem 3.2 is now complete. □ 

Corollary 3.3. Under the conditions of Theorem 3.2 we also have for any k G No 

e -tA g . A 2k + a,k+f (M x ^^ ^ A 2(k+l)+a,(k + l)+| (M x [0)T])> 

e -tA g . A 2k+ Q ,k+f ( M x [ ,T]) -> v^A^+^+f (M X [0,T]) . 
Proof. First note that by definition 

A 2 +Q ,i+f (M x [0)T] ) c ^(A 9jmax ) := {/ G L 2 (M,g) \ A g f G L 2 (M, g)}, 

where A g f G L 2 (M,g) is understood in the distributional sense. The boundary condi- 
tions for the Friedrichs extension A g have been discussed in ([MV11], 2.4). The character- 
ization in ([MV11], Proposition 2.4) implies in particular that any continuous element of 
the maximal domain ^(A gmax ) automatically lies in the Friedrichs domain @(A g ). Hence 
A 2+a,i+f (m x [0,T]) C @(k g ) and similarly 

A 2k+a ' k+ t(M x [0,T]) C 9{A k g ). 

Consequently 

V/ G A 2k+a ' k+ f (M x [0, T]) Vj G N , j < k : A j g e~ tAa * f = e~ tAa * A g f. 

Hence using (3.8) and (3.10) we deduce for / G A 2k+Q ' k+ f (M x [0,T]) and any j < k 

||A,e-' A « * f\\ a = ||e-' A « * A g f\\ a < CVi\\f\\ 2k+a . 

Therefore the statement reduces to k — 1 without loss of generality and we can repeat the 
estimates of Theorem 3.2 ad verbatim taking into account only vector fields X G x~ x V e , so 
that 

/3*(XH) G ( PffPtd )- m -VtTG, 
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where G is a bounded polyhomogeneous distribution on ^£. In particular, comparison 
with (3.2), which encompasses the worst case estimate in case of X = A g , we find that in 
the estimates in Theorem 3.2 we have additional psPtd- We check in local coordinates near 
the front face of ^% that \psPtdG\ < \/i, which provides the additional gain in time in all 
estimates as claimed. □ 



4. The contraction mapping argument 

In this section we outline a contraction mapping argument to prove the existence of local 
solutions to certain quasilinear parabolic equations. For a fixed feasible metric g we study 
the problem 

d t u{t) + A g u = B{u) + Q{u,Vu,Au), u(0) = 0, (4.1) 

where roughly speaking B is a not necessarily linear operator consisting of bounded terms 
of fixed regularity and Q contains the quasilinear terms of lower regularity, which must satisfy 
certain quadratic estimates. We remark that below we have not presented the most general 
possible argument here. Using Schauder estimates for the homogeneous Cauchy problem for 
the heat operator it would be possible to incorporate nonzero initial conditions. We will not 
pursue this here. The main result of this section is 

Theorem 4.1. Consider a fixed feasible metric g and the Cauchy problem (4.1) where the 
operators B and Q satisfy the following mapping properties (fceN) 

(i) B : A 2k +«- k +f (M x [0,T]) -► A 2k+a < k+ f (M x [0,T]) 
(ii) Q : A 2k+a ' k +f (M x [0,T]) -> A2(k-i)+a,(k-i)+f ( M x [ 0;T ]) 

such that for any u £ A 2k+a ' k+ 2"(M x [0,T]) with ||u||2fc+a < P there exists some C(p) > 
such that 



(i) \\Bu - Bv\\ 2 k+a < C(p)\\u - v\\ 2k+a , \\B(u)\\ 2k+a < C ' (p) , 

(ii) \\Qu - Qv\\ 2( k-i) + a < C(p) m&x{\\u\\ 2k+Q , |M| 2fe+Q } \\u - i 



2k+a' 



(iii) ||Q«||2(fc-i)+a < C{p) ||w||5 fc+a - 

Then there is a solution u £ A 2k+a ' k+ 2"(M x [0, T]) to equation (4.1) for some T > 
sufficiently small. 

Proof. We prove the statement using the contraction mapping principle on the Banach space 
A 2k+a ' k+ 2 (M x [0, T]). For positive real parameters p and T to be chosen below, set 

Z^ T := {u £ A 2k+a < k+ f (M x [0,T]) : u(x, 0) = 0, ||«|| 2fc+a < p.} . 

In view of the inhomogeneous heat equation (4.1) we are led to define 

^u := f e - (t ~ s)A « (B(u) + Q{u)) ds. 
Jo 

First observe by the mapping properties of B and Q and from the Schauder estimates in 
the previous section that if u £ A 2k+a ' k+ f (M x [0,T]) then ^/u £ A 2+Q ' 1+ § (M x [0,T]). It 
remains to show ^ is a contraction on Z^. The proof of the contraction property consists 
of two parts. First we must show that \l/ maps Z M) r to itself for some choice of p and T, and 
then we prove the contraction estimate works for the same values of p and T. We begin by 
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assuming that /i < 1. Given u G Z^t we write 



*u = / e" (t - s)As (B(u) + Q{u)) ds 
Jo 

= [ e - (t - s)A « (B(u)) ds + [ e~ 
Jo Jo 

=: (*u)! + (*w) 2 . 



(t-a)A g 



(Q(u))ds 



Applying the Schauder estimate to (\&u)2 and using the properties of Q we find 

\\^u) 2 \\ 2k+a <CK\\u\\l k+a . 
If we take \x < min{ 2^, 1}> then for u G Z^y 

\\^n) 2 \\ 2k+a <CK\\u\\l k+a <CK^ <^. 

This fixes the value of \l. This value of fi persists if we take smaller values of T. Regarding 
(\l/w)i we recall Corollary 3.3. This, combined with the choice /i < 1 implies for u G Z^t 
and T < {fi/2CK) 2 



IK*«; 



1 II 2fe-( 



< ckVt < 



/' 



Combining these two estimates now shows that \P maps the ball Z^t to itself, provided 
/i, T > are sufficiently small. We now check that \l/ is a contraction. For u, v G ^,t 5 we 
must estimate the A 2k+a ' k+ 2"(M x [0,T]) norm of ^u — *£?v. We use the triangle inequality 
to break up the integral into two pieces. Using Corollary 3.3 we find 



f e- {t - s)A °(B(u)-B(v))ds 
Jo 



< CKVfWu 



2k+a 



V \\2k+a< 2 W U - V hk+c 



Similarly, using the quadratic estimates for Q 



-^'^(QM-Qiv^ds 



<CKfi\\u-v\ 



2 AH 



2k+a 



< g \\ U - V hk+a 



We conclude that \1/ is a contraction on Z^t- By the Banach fixed point theorem, \1/ has 



a unique fixed point in Z^t- This concludes the proof of Theorem 4.1. 



□ 



5. Short-time existence of the Yamabe flow 

The Yamabe flow is the following geometric evolution equation 

d t g(t) = -scal(<7(£)) • g(t), g{0) = g im , (5.1) 

where seal (#(£)) denotes the scalar curvature of the metric g(t). See [Brell] for a recent 
survey of results related to this flow on compact manifolds. As an ansatz we will look for 
solutions within the same conformal class as g\ n \ t - In particular, consider 

g (t) = e 2u ^ ■ g mit . 

This transforms the Yamabe flow to a scalar equation for u, where we write A for the 
Laplacian associated to g init 



8 t u = --e- 2u (2(m - l)Au - (m - 2)(m - 1)|Vm| 2 + scal(# init )) , u(p, 0) = 0, (5.2) 
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A few remarks are in order. For a feasible edge metric such as g in \ t it is easy to check that 
in a neighbourhood of dM, the metric takes the form g init = g + h where \h\ go = 0(x 2 ) and 

g = dx 2 + \] ab (y)dy a dy b + x 2 K AB (y, z)dz A dz B . 

The curvature of such a metric is unbounded near the edge x = 0. Straightforward 
computations show that the singular terms in the expansion for scalar curvature of go are 

scal(<7o) = x- 2 (scal(w) - /(/ - 1)) + 0(1). 

Thus in order to make sense of the Yamabe flow we will need to impose further conditions 
on the scalar curvature of the link, namely that scal(ft) = /(/ — 1). This obstruction to the 
Yamabe problem was already noted in [AB03]. With this condition the terms that remains 
are smooth up to the boundary. However, this appears to be insufficient for our contraction 
mapping argument, which requires that various singular derivatives, namely x~ x d z and A So 
of scal(go) be in A a (M,g ), where A a (M,g ) C A a 'f (M x [0,T]) is the subspace of time- 
independent elements. Therefore for the remainder of the argument we assume scal(gi n ; t ) G 
A 2+Q (M, <7i n it) in addition to feasibility. Note that this condition forces scal(«) = /(/ — 1). 

At the time of writing, the authors do not know if there is a more natural geometric 
assumption to impose that would yield the existence of the flow. Fortunately, this assumption 
is satisfied in a variety of interesting cases. For a rigid metric of the form 

g = dx 2 + \) a b(y)dy a dy b + x 2 KAB{z)dz dz B . 

The scalar curvature is exactly 

scai(<7o) = x~ 2 (scal(ft) — /(/ — 1)) + scal(f)). 

Since the metric f) depends on y G B only, it is clear that requiring scal(ft) = /(/ — 1) 
yields scal(g ) G A 2+a (M, g ). Moreover, straightforward but messy computation shows that 
we can also take metrics of the form g = g + h, where g$ is as above and h is a symmetric 
two tensor with \h\ go = 0(x 4 ). To see this compute an expansion for scal(g) keeping track 
of the first two terms in the expansion for g, and observe that requiring \h\ go = 0(x 2 ) forces 
scal(g) to be bounded and smooth with respect to d x , d y , d z down to x — 0. Requiring an 
additional two orders of decay allows this quantity to be bounded in terms of x~ x d x , x^ 1 d z 
and x~ 2 A z . 

We now return to the conditioning of the Yamabe flow equation. Using the exponential 
series we write for any u G A 2k+a,k+ 2"(M x [0, T]) 

e - 2 « = l + n^y4-^=:l + M G( M ), 



n=0 



(n + 1) n\ 



where it is straightforward to check that the elements e~ 2u , G(u) G A 2k+a,k+ 2'(M x [0, T]) are 
locally uniformly bounded. More precisely, for ||«|| 2 fc +Q; < I 1 there exists a constant C(//) > 
such that ||e~ 2 "|| 2fc+a , ||0(u)|| 2fe+Q( < C(/x). Substituting e~ 2u = 1 +uG{u) into the Yamabe 
equation (5.2) we obtain 

d t u + (m - l)Aw = -(m - l)uG(u)Au + ~(m - 2)(m - l)e _2 "|VM| 2 - -e _2u scal(5f init ), 

_ _ 



;i(i 
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d t u + Aw = -uG(u)Au + -(m - 2)e~ 2u \Vu\ 2 - 



-2u 



2(m-l) 



scal(g 



init J i 



Note that by a rescaling of the time variable we may remove the factor of m — 1 multiplying 
the Laplacian. We will instead study 

2 

without further comment. Now set 

1 

t 
1 



Q( u ) ■= - u G{u) Au+-(m- 2)e~ 2u \ Vw| 2 , 



B u) :-- 



-2u 



seal {g 



init ) 



2(m-l) 

In order to apply the short-time existence result of the previous section, we need only 
verify the mapping properties of Q and B. 

Lemma 5.1. Q : A 2k+Q ' k +f (M x [0,T]) -» A2(k-i)+a,(k-i)+f ( M x [o^T]). Furthermore, for 

< ji, there exists a constant C{ji) > 



anyu,u' G A 2k+Q ' k+ f (M x [0,T]) with 
such that 



\ U \\2k+a ' W U \\2k+a 



\\Qu - Qw'|| 2 (fc-i)+a < C(n) max{||M|| 2fc+a , ||w'|| 2fc+a } ||w - u'\ 

||<5w||2(fc-l) < C(fl) ||«||2 fc+a - 



2k+a i 



Proof. It suffices to consider k — 1, the higher order cases treated ad verbatim. The first 
mapping property is clear by the definition of the space A 2+a,1+ 2'(M x [0, T]) and the fact 
that the product of functions in A a '2"(M x [0, T]) is again in A Q '2"(M x [0, T]) by the Leibnitz 



formula. Now if u 6 A 2+a - 1+ 2(M x [0,T]) with 
in A a 'f(M x [0,T]), and 



\u\ 



2+a 



< \x then the first term of Q(u) lies 



(to — l)uG(u)Au\\ a < C(/x) \\u\ 



2+a 



I Ad 



<C(n)\\u\\ 2 2+a . 

For the second term, observe that for u G A 2+a ' 1+ f (M x [0,T]), \Vu\ G A Q 'f (M x [0,T]) 
by the fact that or 1 V c u G A a '?(M x [0,T]). Consequently, 

1 



-(m-2)(m- l)e" 2u |Vw 



< C(n) \\u\ 



2+a 



The second estimate follows from the fact that Q is locally Lipschitz. We estimate only 
one term, the remaining terms are similar. In what follows we omit constants depending only 
on m. First observe that fixing and (p, t) we have 



-2ui 



Vu\ 2 -e 



-2u' 



Vu 



l\2 



-2u 



,-2u'\ 



Vu\ 2 + e~ 2u (\Vu\ 2 - \Vu'\ 2 ) 



-2e~ 2 ^(u - u')\Vu'\ 2 + e- 2M '(|VM| + |V«'|)(|Vu| - |W|), 



where £, lies between u(p,t) and u'(p,t) by the mean value theorem. Estimating this 
expression using the triangle inequality in the A a, 2"(M x [0, T]) norm we find that 



e- 2u \Vu\ 2 -e 



-2u' 



Vu 



/|2 



< C(jj) \\u'\ 

< C'(fi) max{||u|| 



2+a W U U Ib+Q 



2+a ' 



U 



\2+a 



C(/x) max{||w| 

}\\ U - U '\\2+ai 



2+a 1 W U II2- 



\ \\u — u'\ 

-ai 11 



2+a 



as required. 



□ 
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Lemma 5.2. Assume that g init is a feasible edge metric with scal(g init ) G A 2k+a (M,g init ). 
Then B : A 2k+Q ' k+ §(M x [0,T]) -> A 2k+a ' k+ §(M x [0,T]). Furthermore, for any u,u' G 
A 2k+a,k+ 2"(M x [0, T]) with \\u\\ 2k+a , \W\\ 2k+a — ^' ^ere exists a constant C(/x) > such 
that 



\\B(u)-B(u')\\ 2k+a <C(n)\\u-u'\\ 2k+a , 
\\B(u)\\ 2k+a <C(»). 

Proof. The Holder space A 2k+a ' k+ 2'(M x [0, T]) preserved under composition and hence B has 
the desired mapping property. Moreover for any u G A 2k+a ' k+ 2"(M x [0, T]) with HwI^+q < t 1 
we have ||e~ 2n || 2fe+a < C(fj) for some constant C{ji) > 0. This proves the second estimate. 
Regarding the first estimate, this follows as in Lemma 5.1 by the fact B is locally Lipschitz. 

□ 

Using Lemmas 5.1 and 5.2 we may now cite Theorem 4.1 to obtain a short-time solution 
to equation (5.2) and hence equation (5.1). This finishes the proof of Theorem 1.8. 
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